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ABSTRACT 

This paper reports dual-epoch, Very Long Basehne Array observations of H I absorption toward 
3C 147. One of these epochs (2005) represents new observations while one (1998) represents 
the reprocessing of previous observations to obtain higher signal-to-noise results. Significant H I 
opacity and column density variations, both spatially and temporally, are observed with typical 
variations at the level of At « 0.20 and in some cases as large as At « 0.70, corresponding 
to column density fluctuations of order 5 x 10^^ cm~^ for an assumed 50 K spin temperature. 
The typical angular scale is 15 mas; while the distance to the absorbing gas is highly uncertain, 
the equivalent linear scale is likely to be about 10 AU. Approximately 10% of the face of the 
source is covered by these opacity variations, probably implying a volume filling factor for the 
small-scale abso rbing gas of n o more than about 1%. Comparing our results with earlier results 



toward 3C 138 ([Brogan et al.h . we find numerous similarities, and we conclude that small-scale 
absorbing gas is a ubiquitious phenomenon, albeit with a low probability of intercept on any given 
line of sight. Further, we compare the volumes sampled by the line of sight through the Galaxy 
between our two epochs and conclude that, on the basis of the motion of the Sun alone, these 
two volumes are likely to be substantially different. In order to place more significant constraints 
on the various models for the origin of these small-scale structures, more frequent sampling is 
required in any future observations. 

Subject headings: galaxies: individual (3C 147) — ISM: general — ISM: structure — radio lines: ISM 
— techniques: interferometric 



Introduction by iDieter et al. a variety of H I absorp- 



tion studies over the past three decades have found 
Beginning with a two-antenna very long base- AU-scale optical depth variations in the Galactic 

line interferometric (VLBI) observation of 3C 147 interstellar medium (ISM). The initial detections 



1 



were confirmed bv lDiamond et aL ( 19891 ). and the 
first images of the small-scale H I in t he direction 



of 3C 138 and 3C 147 w ere made bv lDavis et al 



( 19961) using MERLI N. iFaison et all (|l998l ) and 



Faison fc Gosd (|200ll) used the Very Long Base- 
line Array (VLBA) to improve the resolution to- 
ward a number of sources to approximately 20 mas 
10 AU). Significant variations were detected in 
the direction of 3C 138 and 3C 147, while no sig- 
nificant variations in H I opacity were found in the 
direction of five other compact radio sources. 

An independent means of probing small-scale 
neutral structures is multi-epoch H I absorp- 
tion meas urements of high proper motion pulsars 



I Fra il et al. 1994 : , Johnston et a l. 2003; Stanimirovic et al 
20031 ). While early pulsar observations suggested 
that small-scale structure might be ubiquitous, 
more recent observations suggest that it could be 
more sporadic. A significant advantage of VLBI 
observations is that they provide 2-D images of 
the opacity variations, rather than 1-D samples as 
in the case of pulsars observations. 

Brogan et al. 1 12005) revisited the observations 



of 3C 13 8, by re-analyzing the 1995 VLBA obser- 
vations ( Faison et al.lll99"8[ ) and by obtaining two 



new epochs of observations (1999 a nd 2002). They 
confirmed the initial results of IFaison fc Gosd 



(|200l[ ). that there are small-scale opacity changes 
along the line of sight to 3C 138 at the level of 
Axmax = 0.50 ± 0.05, with typical sizes of roughly 
50 mas 25 AU). However, with multiple epochs 
and improvements in data analysis techniques 
(yielding an increase of a factor of 5 in the sensi- 
tivity of the 1995 epoch), they reached a number 
of additional significant conclusions: 

1 . They found clear evidence for temporal vari- 
ations in the H I opacity over the seven- 
year time span of the three epochs, consis- 
tent with structures moving across the line of 
sight at velocities of a few tens of kilometers 
per second, though the infrequent sampling 
in time means that they could not determine 
whether these structures were persistent. 

2. They found no evidence for a drop in the H I 
spin temperature, as would be evidenced by 
a narrowing of line widths at small scales 
compared to single dish measurements. In 
turn, a constant H I spin temperature im- 
plies that the small-scale opacity variations 



are due to density enhancements, although 
these enhancements would necessarily be ex- 
tremely over-pressured relative to the mean 
interstellar pressure, far from equilibrium, 
and likely of relatively short duration. 

For the first time they determined that the 
plane of sky covering fraction of the small- 
scale H I gas is roughly 10%. In turn, this 
small covering fraction suggests that the vol- 
ume filling factor of such gas, within the cold 
neutral medium, is quite low (< 1%), in 
agreement with HST obse rvations of high- 
pressure gas in th e ISM (|Jenkins fc Tripp 
200ll:|Jenki^l2004 . 



4. They simulated pulsar observations that 
have been used to search for H I opacity 
variations and showed that the existing pul- 
sar observations have generally been too 
sparsely sampled (in time) to be useful in 
studying the details of small-scale H I opac- 
ity variations. 



Wh ile the multi-epoch study of iBrogan et al 
(|2005l ) represented a substantial improvement 



nonetheless their conclusions rested on observa- 
tions of only one line of sight. In light of this 
sample of one, their conclusions might seem rather 
audacious , particularly given the la rger sample ob- 
served by IFaison et al. I (|l998h and IFaison fc Gosd 
( 200ll ). in which most of the objects did not show 
variations in the H I absorption. The 3C 138 
study has shown that the key to a successful small 
scale H I study is a background source with both 
high surface brightness (> 60 mJy beam~^) and 
large angular extent (> 100 mas). The quasar 
3C 147 is one of the few sources that shares these 
characteristics with 3C 138. This paper presents 
dual-epoch observations of 3C 147 that were de- 
signed specifi c ally t o confront the conclusions of 
Brogan et al. ( 20051 ) with a second hue of sight. 
Section [2] of this paper describes the observations, 
focussing on the new observations acquired for 
the second epoch, discusses the results, and 21 
presents our conclusions and recommendations for 
future work. 

2. Observations 

We have observed the Galactic H I absorp- 
tion (near —10 km s~^) toward the quasar 3C 147 
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at two epochs. Epoch I was 1998 October 22, 
and the resuhs from those observations ha ve been 
pubhshed previously by iFaison fc Gosd ([2001). 
Epoch II consists of new data observed on 2005 
August 21. Table [T] summarizes the basic observ- 
ing parameters for the two epochs. 

For both epochs the data were obtained using 
the 10 antennas of the Very Long Baseline Array 
combined with the Very Large Array with its 27 
antennas operating in a phased-array mode. For 
the 2005 epoch, the Green Bank Telescope was 
also used. The observing duration was 12 hr for 
the 1998 epoch and 16 hr for the 2005 epoch, in- 
cluding time spent on calibration sources. The 
proximity of the VLA to the VLBA antenna at 
Pie Town, New Mexico, significantly increased our 
sensitivity to large-scale structures. Four sepa- 
rate spectral windows or intermediate frequency 
bands (IFs) were used, with one IF centered on 
the absorption line (at an approximate LSR veloc- 
ity of —10 km s~^) and three IFs separated by at 
least 100 km s^^ in velocity in order to sample the 
21 cm continuum emission. For the 1998 epoch, 
the data were correlated with velocity channels 
of 0.4 km s~^ , with a bandwidth of 500 kHz per IF 
over 256 spectral channels; for the 2005 epoch, im- 
provements in the correlator allowed the number 
of spectral channels per IF to be increased to 512, 
with a concomitant improvement in the velocity 
resolution to 0.2 km s~^. 

Broadly similar data reduction procedures were 
used for the two epochs. For the 2005 epoch, 
the data were calibrated for the frequency depen- 
dence of the bandpass using observations of 3C 48 
and amplitude calibrated using system tempera- 
tures measured at the individual antennas. The 
most significant difference in the calibration is 
that for the 2005 epoch, we attempted to phase- 
reference the observations to the compact source 
IVS B0532+506, separated by 1?3 from 3C 147. 

Our initial motivation for this change in pro- 
cedure is that 3C 147 has a sufficiently com- 
plex structure that fringe-fitting assuming a point- 
source model could yield erroneous residual delay 
and rate solutions. In practice, phase referencing 
did not prove useful. The phase-referencing cycle 
time was short enough that latency in the VLA 
system often resulted in the VLA acquiring no 
data. The most significant difficulty, however, was 
that only one of the epochs was phase-referenced. 



There was an apparent ofi^set in the core position 
between the two epochs (with a magnitude of a 
fraction of the synthesized beam width or a few 
milliarcseconds) that biased any attempt to com- 
pare results from the two epochs (e.g., compar- 
ing the integrated line profiles). Consequently, we 
did not make use of the phase-referenced data for 
constructing the H I line profile or opacity im- 
ages. One obvious impact on our results is that 
the sensitivity of the 2005 epoch H I line data 
is less than it could have otherwise been due to 
the phase-referencing cycling between 3C 147 and 
IVS B0532-h506. 

Two of the three continuum IFs were then av- 
eraged together and several iterations of hybrid 
imaging (iterative imaging and self-calibration) 
were performed. After the final iteration of self- 
calibration, the phase and amplitude solutions 
were applied to the IF containing the H I line. 
The line-free velocity channels in this IF were av- 
eraged together to produce a continuum data set, 
which underwent a final round of hybrid imaging, 
the solutions from which were applied to the ve- 
locity channels containing the line. Finally, the 
continuum emission was subtracted from the ve- 
locity channels containing the H I line and the 
resulting line data set was imaged. 

We also reprocessed the observations of Faison fc GossI 
(|200ll ) in a similar fashion. A si gnificant differ- 
ence f rom the original analysis of iFaison fc Gosj 
(|200l[ ) is that we used the continuum image from 
the 2005 epoch as an initial model for fringe fitting 
the 1998 epoch data (for which no phase referenc- 
ing was performed). The combination of a better 
initial model and improvements in the imaging 
software and analysis procedures led to a substan- 
tial improvement in the reprocessed 1998 epoch 
data. The noise in the 1998 epoch continuum im- 
age has improved by nearly an order of magnitude, 
and the improvement in the spectral line images 
is a factor of a few. As was the case for 3C 138, 
the original analysis found a significantly higher 
peak brightnes s than we do , by a similar factor 
(w 30%). Like iBrogan et all (|2005h . we attribute 
thi s difference to t he us e of a point source model 
bv IFaison fc GossI (l200lh in the original fringe fit- 
ting along with other details of the subsequent 
imaging and self-calibration. 

Faison fc Gosd 



Fol lowing the procedure of 



(|2001l ). both the continuum images and continuum- 
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Table 1 
Observational Log 



Parameter 


Faison & Goss ("2001") 


Epoch I 


Epoch II 


Date 


1998 October 22 


1998 October 22 


2005 August 21 


Number of IFs 


4 


4 


4 


Bandwidth per IF (MHz) 


0.5 


0.5 


0.5 


Spectral channels 


256 


256 


512 


Channel separation (km s~^) 


0.41 


0.41 


0.21 


Velocity resolution (km s"^)" 


0.4 


0.41 


0.21 


CLEAN beam (mas)'' 


5x4 


8.2 X 5.6 


7.6 X 7.1 


Continuum Peak (Jy beam"^)'^ 


2.31 


1.681 


1.801 


Continuum rms noise (mjy beam"^)'^ 


6.5 


1.7 


1.1 


Spectral line rms noise (mJy beam"^)'^ 


7.6 


5.0 


3.5 


General Parameters for 3C 147 


Position (equatorial, J2000) 


05'^42™36?13788 


+49°07'51'.'2335 




Position (Galactic, longitude & latitude) 


+161.69° 


+10.30° 




Redshift 


0.545 







''After Hanning smoothing during imaging process. 

''The CLEAN beam before convolution to 10 mas. All subsequent values are for the convolved 10 mas 
resolution images. 



'^For the image after it has been convolved to 10 mas resolution. 



Note. — The values listed under the lFaison &: Gos3 ll200lh column are for the original analysis. The values 
listed under the Epoch I column are for this analysis, after the reprocessing of the data as described in the 
text. 
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subtracted line cubes were convolved to 10 mas 
resolution. The u-v coverages for the visibility 
data from the two epochs were similar, producing 
images with angular resolutions of approximately 
7 mas (Table [T]) . The convolution of the contin- 
uum images and continuum-subtracted line cubes 
is an attempt to minimize any effects of modest 
differences in the u-v coverage between the epochs. 
The second epoch data were also smoothed in ve- 
locity so that their velocity resolution matched 
that of the first epoch. 

An optical depth cube, calculated as th i{a, S, v) 
-ln[l - /iino(a,(5,f)//cont(a,<5)], where hnc and 
-^cont are, respectively, the images from formed 
from line and line-free channels. Because the 
signal-to-noise ratio in the optical depth images 
is low where the continuum emission is weak, the 
optical depth images were blanked where the con- 
tinuum emission was less than 5% of the peak 
emission. 

3. Results 

3.1. 21 cm Continuum 

Figure [T] presents the 21 cm continuum im- 
age of 3C 147 from the new observations of the 
2005 epoch. There is good qualitative agree- 
ment between our image and previously published 
images at comparable wave l engths (18-20 cm , 



Readhead fc \yilkinsonl 1980l: Zhang et a 



1991 



Polatidis et al.l 119951 : iFaison fc Gosd I2001I ) . The 

source displays its well-known core-jet structure, 
with the jet extending some 200 mas to the south- 
west before bending to the north. Also promi- 
nent is diffuse emission to the east of the core, ex- 
tendin g to the north, first noticed bv lZhang et al. 
(|l99l[ ). The resolution of our observations is not 



high enough to resolve t he northeast extension 
from the core found by iReadhead Sz Wilkinson 
(|l98flf ). 

We assessed the continuum images from the two 
epochs for variability. Creating a difference image 
between the two epochs, we find that any vari- 
ability in the source is below the 20 mJy beam~^ 
level. Even if the source is variable at this level, 



as for the iBrogan et al.l (j2005i ) analysis, variabil- 
ity will not impact our optical depth calculations, 
because (1) the continuum appropriate for each 
epoch was used in the optical depth calculations, 
(2) amplitude self-calibration solutions were never 





-100 -200 

MilliARCSEC 



Fig. 1.— The Epoch II (2005) 21 cm contin- 
uum image of 3C 147 obtained with the VLBA 
and phased VLA. The clean beam is 7.6 mas 
X 7.1 mas at a position angle of —15°. The 
rms noise level is 0.5 mJy beam^^, and the 



contours are set at 0.5 mJy beam 



-3, 5, 



7.07, 10, 14.1, 20, . . .. The gray scale is lin- 
ear between 1.5 and 500 mJy beam"-'^. This 
image shows the source at the full resolution; 
for subsequent analysis, the image was convolved 
to 10 mas resolution. The origin is at (J2000) 
OS'' 42'° 36!1379 +49° 51' 07'.'234. 



transferred between the epochs, and (3) the intrin- 
sic continuum morphology at 10 mas resolution 
does not appear to change from epoch to epoch. 

The flux density in our image is 18 Jy. The VLA 
Calibrator Manual lists of flux density of 22.5 Jy, 
indicating that our observations have recovered 
80% of the source's total flux density. 

3.2. H I Line Profile 

Figure [2] shows the average optical depth pro- 
flle from the new observation of 2005. Even 
with its relatively high latitude (6 = -flO°), 
the prof ile is comple x, making it similar to 



3C 138 (|Brogan et al.ll2005l) . Our profile is in 



good agreement with previously published p rofiles 
(|Kalberla etld] Il985l: [Faison fc Gos3 l200ll ). and 
the difference between the line profiles from the 
two epochs sh o ws on ly modest variations. Like 
Faison fc Gosd ( 200ll ). we shall restrict our atten- 
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-20 20 40 

Velocity (km/s) 

Fig. 2. — The solid line shows the average H I op- 
tical depth profile toward 3C 147 from Epoch II 
(2005). The dotted line shows the difference be- 
tween the line profiles from the two epochs. Also 
marked are the three significant velocity compo- 
nents at which further optical depth analysis is 
performed (viz. Tabled]). 



tion to the three most prominent velocity compo- 
nents, those with approximate central velocities 
of -10.4, -8.0, and 0.4 km s^^ (see below). 

Faison fc Gosd (2001) have discussed the dif- 



ficulties with assessing the distance to the ab- 
sorbing gas. Under the simple assumption that 
all H I gas is confined to a 100 pc thick layer, 
the absorbing gas must be within 0.6 kpc. How- 
ever, the kinematic distance to the gas causing 
the —8 km s^^ absorption is uncertain, with dis- 
tances as large as 1.1 kpc allowed. As a nominal 
value, we adopt the conversion that our resolu- 
tion of 10 mas corresponds to a linear distance 
of 7.5 AU (1 mas = 0.75 AU), implying a distance 
of 750 pc to the gas, though differences of as much 
50% are possible. 

For subsequent analysis, we fit the 2005 epoch 
optical depth line cube by gaussian components, 
using the profile of Figure [2] as an guide to ini- 
tial values for the component parameters. For 
the fitting, we focussed on the three significant 
components identified. The fitting was done on 
a pixel-by-pixel basis, with an independent three- 
component fit for each pixel. Table [2] summarizes 



results of the fits, averaged over the face of the the 
source. Guided by the results of the fitting from 
the 2005 epoch, a similar fitting was performed for 
the 1998 epoch. 

3.3. Small-Scale Structure 

Figures [3] and |4] show column density fluctua- 
tion images at the two epochs for the three differ- 
ent velocity components. From the gaussian fits, 
column density images were constructed from the 
fitting results by multiplying the maximum optical 
depth by the velocity width (iVni/Ts oc T(T„, see 
below regarding the spin temperature T^). In or- 
der to highlight fluctuations, the average column 
density from the 2005 epoch across the face of the 
source was subtracted from these column density 
images to produce the column density fluctuation 
images. The signal-to-noise is not uniform across 
the face of the source, and tends to decrease near 
the edges. In order to aid in assessing the reality 
of features. Figures [3] and |4] also show the column 
density fluctuation signal-to-noise ratio images. 

We show column densi t y fluc tuation images, in 
contrast to iBrogan et al. ( 20051) who showed op- 
tical depth channel images. For 3C 147, analysis 
of the column density fluctuations is required be- 
cause the velocity fleld of the absorping gas ap- 
pears to change some what within each ve locity 
component. For 3C 138,|Brogan et al.l(|20n5h com- 
pared the optical depth channel images at differ- 
ent velocities and concluded that any velocity field 
fluctuations were negligible, in contrast to the sit- 
uation for 3C 147. Also, in converting to column 
density, we assume a uniform spin temperature of 
the gas of = 50 K (|Heileslll997l) . Unlike 3C 138 
which was part of the Millennium Arecibo 21 cm 
Survey (|Heiles fc Trol and 2003), 3C 147 is outside 
the Arecibo declination range so that there is less 
direct information about the absorbing gas along 
this line of sight. 

Figures [IHS] show cuts through the column den- 
sity fluctuation images. In all cases, spatially sig- 
nificant changes in the column density between 
the two epochs are clearly apparent for all of the 
H I components, at significance levels exceeding 5ct 
over most of the face of the source. We have also 
conducted an analysis in which we consider a con- 
stant column density cross-cut to be the null hy- 
pothesis. In a sense, the null hypothesis can be 
firmly rejected as typical values, for both epochs 



6 
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.100 -120 -140 -160 -180 



Fig. 3. — Column density fluctuation images, for the 1998 epoch (Epoch I), with the average value of the 
2005 epoch column density subtracted, ANm/Tg oc J dv{T — (r2oo5))- Horizontal white lines are pixels 
where the gaussian fitting failed to converge. Left panels show the signal-to-noisc ratio, on a linear scale 
with the dynamic range restricted to 0-15 and with a single contour showing a signal-to-noise ratio of 5. 
Right panels show the column density fluctuations, on a linear scale. The gray scale bar shows the column 
density fluctuation, in units of 10"'^^ cm~^. (Top) Column density fluctuations for —10.4km s"-'^ velocity 
component in 1998, the 2005 average optical depth that was subtracted is (T2005) = 0.31; (Middle) Column 
density fluctuations for — S.Okms"^ velocity component in 1998, (T2005) = 0.93; and {Bottom) Column 
density fluctuations for 0.4 km velocity component in 1998, (T2005) = 0.54. 
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Table 2 

Optical Depth Profile Gaussian Component Fit Results 



Component 


Epoch 


Central Velocity 


Velocity Width 


Maximum Optical Depth 






(km s-i) 


(km s^^) 




1 


1998 


0.4 


5.1 


0.5 




2005 


0.3 


4.9 


0.5 


2 


1998 


-8.0 


1.6 


0.8 




2005 


-8.0 


1.5 


0.8 


3 


1998 


-10.4 


6.1 


0.3 




2005 


-10.4 


6.2 


0.3 



" i 






cm s 


2005 - 












" i 








1998 






1 TiiiH.BH.'^ 






" 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 


1 1 1 1 1 




1 1 " 






i -8.0 


krr, 


2005 - 










i 












... ' Vv 


a..*...i rf 




1998 : 


i 

- 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 


'..'ig,EE!aai 

1 1 1 1 1 




i 




,SHi, -10.4 


km s~' 


1 1 : 










2005 J 








1998 ■ 



50 100 150 200 250 

Distance along cross-cut (mos) 

Fig. 7. — Comparison of a cross-cut along the ma- 
jor axis at the two epochs for the three velocity 
components. For clarity, the 1998 epoch is shifted 
(by 15 X 10^^ cm~^) relative to the 2005 epoch. 
For both epochs, a spin temperature T, = 50 K is 
assumed. Also shown are uncertainties (±2cr), al- 
though in many cases they are only slightly larger 
than the symbol size. Further, because the restor- 
ing beam can induce correlations, we plot only ev- 
ery fifth datum. The major axis illustrated is the 
central of the three in Figure \5\ 

and all velocity components, are ^ 10 (reduced 
X^)- Typical angular scales for column density 
variations are approximately 15 mas, correspond- 
ing to a linear scale of approximately 10 AU. 
Brogan et al. ( 20051 ) discussed the possible sys- 



tematic effects that might affect the extraction of 
reliable optical depth or column density variations 
from dual-epoch imaging such as presented here. 
We do not repeat their discussion, but consider 



l4 km s 



b8' 



A 



.iiBiii 



mmi 
f- 



— 1 h- 

i.O km s 



jjliHfiii 



10.4 km s ^ 



,iii'=.!^,iSiHHiaa.ii.' ''ii' % 



20 40 60 80 100 

Distance along cross-cut (mas) 

Fig. 8. — As for Figured (and Figure [6]), but for 
a minor axis slice. 

many of the same issues and conclusions to hold. 
Namely, while small differences in the images from 
epoch to epoch may be due to the details of the 
observations and data reduction (e.g., spatial fre- 
quency or u-v coverage, slight differences in the 
imaging), a number of steps were taken during 
the analysis in an effort to minimize the differ- 
ences between the epochs, and the column density 
variations are significant. 

The column density fluctuations in Figures [SI- 
IT] correpond to peak-to-peak optical depth vari- 
ations as large as At w 0.7, and typical opti- 
cal depth variations on scales of approximately 
15 mas ranges from 0.1-0.3. The associated un- 
certainties in the optical depth are a-r ~ 0.07, 
implying significant variations at the 3a level, 
and the column density cross-cuts indicate varia- 
tions at even higher significance are present. Fur- 
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MilliARC SEC 



Fig. 5. — {Top Left) Cross-cuts taken approximately along the major axis showing the column density 
fluctuations for the 1998 epoch. ( Top Right) Cross-cuts taken approximately along the major axis showing 
the column density fluctuations for the 2005 epoch. (Bottom) Illustration showing where the cross cuts were 
taken taken, with the column density fluctuations from the 0.4 km s~^, 2005 epoch shown for reference. 
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ther, the magnitude of the variations is approxi- 
mately correlated with the strength of the aver- 
age H I absorption toward 3C 147, in that the 
largest variations are observed at —8.0 km s~^, 
followed by 0.4 km s~^, with the smallest varia- 
tions at -10.4 km s^^ (cf. Figures [2] and [5]-[8l) . 

The opacity variations toward 3C 138 show 
changes that are consistent w ith motion of struc- 
tures across the line of sight (|Brogan et al.ll2005i) . 
though there is considerable uncertainty with 
making these identifications (as they discuss). 
Possible motions of structures across the line 
of sight toward 3C 147 are also visible in the 
cross-cuts. Examples of such motions include the 
features at distances between 50 and 100 mas 
(Figures [7]). Caution in interpreting these fea- 
tures as arising from motions is clearly warranted, 
however, given that we have only two epochs. 
Nonetheless, typical position shifts appear to be of 
the order of 5 mas. Over the 7 yr interval between 
observations, the implied proper motion is just un- 
der 1 mas yr^^, equivalent to a velocity of order 
3 km s~^, at a distance of 750 pc. For comparison, 
and recognizing that there i s considerable uncer- 
tainty in these comparisons, Brogan et al. ( 20051 ) 
find larger values for the apparent velocities of 
structures toward 3C 138 (« 20 km s~^). 

One of the motivations for undertaking these 
observations was to assess whether the optical 
de pth variations, both s patial and temporal, found 
by iBrogan et al. toward 3C 138 indicated 

that the line of sight to that source was in some 
sense "special" or anomalous. Comparison of Fig- 
ures [3] and [4] and Figures [SHHl with the corre- 
sponding ones from lBrogan et all ( 2005h show that 
they are qualitatively similar, with clearly signif- 
icant opacity or column density variations occur- 
ring both in space and time. 

Quantitatively, there are modest differences in 
the opacity/column density variations between 
3C 147 and 3C 138. We estimate that the typi- 
cal angular scale of opacity variations is 15 mas 
(« 10 AU), as opposed to about 50 mas (~ 25 AU) 
toward 3C 138, though the linear scales are com- 
parable. The magnitude of the variations toward 
3C 147 also seems somewhat smaller. Optical 
depth changes (both in space and time), and cor- 
responding column density fluctuations, are a fac- 
tor of a few to several larger for 3C 138 — optical 
depth changes of 0.4 and larger (> 10^" cm~^) for 




2005. -10.4 km s 



2005, -a.O km s 



2005. 0.4 km s 



At = |t - <T>| 



At = |t - <T>| 



Fig. 9. — Fractional number of pixels in an optical 
depth channel image as a function of the optical 
depth variation At = |t — (r)|. Left panels show 
the first epoch (1998), and right panels show the 
second epoch (2005). (Top) -10.4km s'^; (Mid- 
dle) -8.0 km s"^ and (Bottom) 0.4 km s"^ 



3C 138 vs. optical depths typically not exceeding 
0.3 (- 5 X 10^9 cm"2) foj. 3Q i^-j 

3.4. Small-Scale H I Covering and Filling 
Factors 



Brogan et al.l (|2005l ) used their observations of 
3C 138 to conclude that the (two-dimensional) 
covering factor of small-scale H I opacity varia- 
tions was about 10%, from which they inferred 
a three-dimensional filling factor of probably less 
than 1%. Although the optical depth variations 
appear qualitatively similar for lines of sight to- 
ward 3C 138 and 3C 147, we have repeated their 
analysis in order to determine the covering and 
filling factors for the line of sight to 3C 147. 

Figure [9] shows the fractional number of pixels 
in a optical depth channel image for the three dif- 
ferent velocities at both epochs. Most of the opac- 
ity variations are at a level less than about 0.2 in 
optical depth, and we do not consider them signifi- 
cant. Restricting to optical depth variations larger 
than approximately 0.2 (w 3cr), most of the cov- 
ering fractions are about 10%, ranging from a low 
value of a few percent to a high exceeding 25%. 

While we have decomposed the optical depth 
profile into gaussian components, our decomposi- 
tion is likely not unique and we fit only for three 
components, whereas Figure [2] clearly shows that 
there are could be more components. Thus, a 
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plausible upper limit to the volume filling factor 
of the small-scale absorbing gas is obtained by as- 
suming that there are multiple components, each 
of which contributes equally. We obtain aii upper 
limit of 1%, a value which, as Brogan et al. ( 20051 ) 
emphasize, is not directly measurable. 

4. Discussion and Conclusions 

We have presented two epochs of observations 
of the H I optical depth across the face of the 
source 3C 147 on scales of approximately 10 mas. 
The motivation for these observations was assess- 
ing whether spatial and temporal H I opacity vari- 
ati ons found in n i ulti-ep och observations of 3C 138 
bv iBrogan et al. ( 20051 ) were in some sense "spe- 
cial" or anomalous. 

We find qualitatively similar opacity and col- 
umn density variations toward 3C 147 as were 
found toward 3C 138. Quantitatively, the vari- 
ations toward 3C 147 appear to be somewhat 
smaller in angular scale (15 mas vs. 50 mas) and 
smaller in magnitude (by a factor ~ 5). While 
the typical angular scale of the variations toward 
3C 147 appears smaller, the absorbing gas may be 
more distant than that causing the absorption to- 
ward 3C 138 If so, the resulting linear scales 
are comparable (10 AU for 3C 147 vs. 25 AU for 
3C 138), though the uncertainties are large. 

Further similarities are observed in the covering 
and filling factors of the small-scale absorbing gas 
toward both sources. For both lines of sight, the 
covering factor appears to be approximately 10%, 
and the volume filling factor, while not measured 
directly, has a plausible upper limit of 1% (and 
potentially much less). 

Both 3C 138 and 3C 147 display significant 
H I opacity variations across their faces, im- 
plying variations within the ISM on scales of 
about 10 to 50 AU over path lengt hs ranging 



from 100 to 1000 pc. Using MERLIN, Goss et al 



(|2008l ) also have resolved significant H I opac- 
ity variations across the faces of 3C 111, 3C 123, 
and 3C 161, implying structure on scales of 50 
to 500 AU. We conclude that the conditions 
that cause such small-scale variations are fairly 
widespread within the Galactic ISM. The reason 
that so few other sources display su ch small-scale 



and 3C 147 have the combination of angular ex- 
tent and surface brightness required to conduct 
these observations. One unfortunate implication 
of this conclusion is that milliarcsecond-scale H I 
observations of other sources will largely not be 
useful for probing the small-scale structure with- 
out a significant increase in sensitivity. One pos- 
sible target, particularly with the existing High 
Sensitivity Array (HSA)0 may be 3C 380. 

In one aspect, however, the lines of sight to 
3C 147 and 3C 138 d o diffe r. For the 3C 138 
analysis, iBrogan et al. ( 2005 ) used optical depth 
velocity channel images whereas, for 3C 147, we 
fit the optical depth line cube with gaussian com- 
ponents and used the resulting column density im- 
ages. This difference in approach was motivated 
by the velocity structure that was apparent within 
an H I component in the 3C 147 optical depth line 
cube. 

We have not been able to find a ready expla- 
nation for this difference. Both sources are seen 
toward the Galactic anticenter, with Galactic co- 
ordinates (longitude, latitude) of (161?7, 10?3) for 
3C 147 and (187?4, -11?3) for 3C 138, respec- 
tively. To first order, the lines of sight to both cut 
almost perpendicular to the Perseus spiral arm. 
Further, were velocity crowding the explanation, 
it would seem that that should be more of an issue 
for 3C 138 than for 3C 147. We have also consulted 
the WHAM Ha survey a nd the Green supernova 
remnant (SNR) catalog (Green 20061 ). reasoning 
that Ha and SNRs might serve as a tracers of 
turbulence injected by winds or explosions from 
massive stars. There are no obvious indications 
that the line of sight to 3C 147 should be affected 
any such turbulence — indeed a comparison of the 
3C 138 and 3C 147 lines of sight suggest that the 
line of sight to 3C 138 would be more likely to be 
the one that would display any such evidence of 
turbulence. 

An alternate possibility is that the kinematic 
differences between these two lines of sight re- 
fiect small-scale fe atures, and possibly the past 
history of the gas. iKalberla et al. imaged 
the H I emission around the line of sight toward 
3C 147 at 1' resolution (« 1 pc linear scale). They 



opacit y variations is likely to be, as IBrogan et al 



(|2005[ ) discuss, that few sources other than 3C 138 



^ The VLBA combined with other large aperture tele- 
scopes such as the phased VLA, the Green Bank Telescope, 
Arecibo, or the 100-m Efflesberg telescope. 
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found a series of filaments and small clumps of H I 
emission, and they were able to associate at least 
some of the absorption features with small emis- 
sion clumps. The amount of small-scale structure 
(in H I emission) toward 3C 147 is not generally 
observed in the o n-going GALFA H I survey at 
Arecibo. Further. [Kalberla et alJ ( 1985 ) find that 
a large fraction ('^ 80%) of the H I in emission in 
the direction of 3C 147 has a temperature of 500- 
2000 K. At this temperature, the gas would be 
thermally unstable. While the warm H I is not 
responsible for the absorption, the possibility that 
this line of sight contains thermally unstable H I 
is consistent with a scenario in which the micro- 
physics, and potentially the past history of the 
gas, leads to kinematic variations within an H I 
absorption component. 



Kalberla et al.l (|l985l ) also suggested a relative 
distance ordering of the gas. Comparison of the 
opacity variations (Figures EHSl) suggests that the 
opacity variations are smallest in amplitude for the 

— 10.4 km s^^ velocity component and increase in 
magnitude for the 0.4 km s~^ and the —8.0 km s~^ 
velocity components. One interpretation is that 
the opacity variations result from structures of 
essentially constant size, which are comparable 
to or smaller than the equivalent linear size of 
our beam (~ 10 AU). If this were the case, we 
could obtain a relative distance ordering of the 
gas, with the —8.0 km s~^ material being the near- 
est, followed by the 0.4 km s^^ material, and the 

— 10.4 km s~^ material being the most distant. 
An alternate interpretation (see above) would at- 
tribute these differences to the history of exposure 
of the gas to shocks or other interstellar distur- 
bances. Whichever is the case, there is clearly 
significant structure on large scales (~ 1 pc), sug- 
gesting that such structure could persist to smaller 
scales. The combination of VLBA and VLA data 
(as well as potentially MERLIN data) might be 
able to explore the connection between the small- 
and large-scale opacity variations. 

iBregman et al. (1983) have set an upper limit 
(3cr) on the magnetic field toward 3C 147 of 
< 50fiG, based on Zeeman effect measure- 
ments in H I spectra. Under the standard assump- 
tion that discret e H I structu res require densities 
n ~ 10 cm~ 3 (lHeiledll997l) . with a typical ve- 
locity width of w « 3 km s~^ (Figure [2] and Ta- 
ble [21), one concludes that magnetic and turbulent 



equipartition requires a magnetic field strength of 
order 400 /xG, well above the observed upper limit. 
As for the optical de pth variations toward 3C 138 
( Brogan et al. the optical depth variations 

cannot be in magnetic and turbulent equilibrium, 
unless there is significant blending and dilution of 
the magnetic field on the angular scales over which 
the Zeeman effect measurements were made. 

Our observations of 3C 147 do not produce 
any new constraints on the nature of the small- 
scale structures vis-a-vis whether they repre- 



sent "statistical" fluctuations (e.g., iDeshpande 



20001) ■ "non-equilibrium" physical structures (e.g . , 
Jenkins k Trippll200ll : [Hennebelle fc Auditll2007l) 



or di screte "tiny scale atomic structures" ( HeilesI 
19971 ). While the level of opacity variations to- 
ward 3C 147 are lower than those toward 3C 138, 
we believe that this lower level can be accommo- 
dated easily within any of these scenarios. A wide 
range of opacity variations might be expected if 
these result from non-equilibrium processes, par- 
ticularly because the level of opacity variations 
co uld depend up on the history of the gas. Also, 
as lBrogan et al. I (12005) note, the predicted level of 
opacity variations within a statistical description 
depends sensitively upon the assumed spectral 
index of the underlying power law; within the cur- 
rent uncertainties for this spectral index, a large 
range of opacity variations is allowed. 

What would be required in order to place sig- 
nificant constraints on these small-scale opac- 
ity variations? Ideally, one should monitor the 
same volume of gas and determine how the struc- 
tures evolve. In the most simple comparison, dis- 
crete structures should show only linear motion 
across the line of sight, while fluctuations or non- 
equilibrium physical conditions might also cause 
the appearance of the opacity variations to change 
significantly. 

The elapsed times betw een observatio n s pre- 



sented here and those in iBrogan et al. (|2005l) 



range from 3 to 7 yr. These intervals are only a few 
percent of the estimated time for discrete struc- 
tures to change substantially (« 500 yr, HeilesI 
20071 ). However, on milliarcsecond scales, the 
line of sight to 3C 147 effe ctively samples a vol- 
ume through the G alaxy (jMarscher et al ] I1993I: 
Dieter-Conkiinll2009l ). The Sun's velocity through 
space causes this volume to move between our 
observing epochs (Figure [TO]) , in addition to any 
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Fig. 10. — An illustration of the different volumes 
of the Galaxy sampled by multi-epoch VLBI ob- 
servations. The distance to the absorbing gas is 
assumed to be 750 pc, the space velocity of the 
Sun is assumed to be 30 km s~^, and the gas is 
assumed to be stationary. The resulting effective 
proper motion is 8.4 mas yr~^; a smaller assumed 
distance for the gas would result in a larger proper 
motion while a smaller space velocity for the Sun 
would result in a smaller proper motion. Shown 
is the apparent position of the source at three hy- 
pothetical epochs, each separated by 3 yr, compa- 
rable to the typical separation in epochs for the 
existing multi-epoch VLBI H I absorption obser- 
vations. (Epoch I is red, Epoch II is green, and 
Epoch III is blue.) The white areas indicate the 
only sampled volumes of gas common to all three 
epochs. 



terval is no longer than about 9 months, with even 
more rapid sampling desirable. 
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tional Science Foundation operated under coop- 
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Basic research in radio astronomy at the NRL is 
supported by 6.1 NRL Base funding. 



motion that the gas itself might have. A sim- 
ple estimate of the Sun's motion suggests that 
an entirely new volume through the Galaxy could 
be sampled on time scales of approximately 3 yr. 
That is, in the typical interval between VLBI ob- 
servations, essentially an entirely new volume of 
the Galaxy is sampled by the line of sight. 

Consequently, we also conclude that the time 
sampling of the existing multi-epoch VLBI obser- 
vations has been too coarse to distinguish between 
the various models for the small-scale opacity vari- 
ations. Ideally, one would like to monitor the same 
volume of gas, to determine if the opacity varia- 
tions appear to be simply in motion or also chang- 
ing in appearance. We estimate that, for either 
3C 138 and 3C 147, an appropriate sampling in- 
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